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Abstract 

We propose a solution of the naturalness problem in the context of the mul- 
tiverse wavefunction without the anthropic argument. If we include microscopic 
wormhole configurations in the path integral, the wave function becomes a su- 
perposition of universes with various values of the coupling constants such as 
the cosmological constant, the parameters in the Higgs potential, and so on. We 
analyze the quantum state of the multiverse, and evaluate the density matrix of 
one universe. We show that the coupling constants induced by the wormholes 
are fixed in such a way that the density matrix is maximized. In particular, 
the cosmological constant, which is in general time-dependent, is chosen such 
that it takes an extremely small value in the far future. We also discuss the 
gauge hierarchy problem and the strong CP problem in this context. Our study 
predicts that the Higgs mass is nih = 140 it 20 GeV and = 0. 
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1 Introduction and Conclusion 



One of the major problems of particle physics and cosmology is the smallness of the 
observed value of the vacuum energy, that is the cosmological constant A. We must 
explain why A is many orders of magnitude smaller than the Planck scale [1]. One of 
the most promising attempts to solve this problem is the one based on the Euclidean 
wormhole effect first proposed by Coleman pjfand studied closely by other authors 
[IHTl]. In this paper, we discuss the wormhole effect in the context of the Lorentzian 
multiversqj, and propose a mechanism to solve the naturalness problems such as the 
cosmological constant, the gauge hierarchy, and the strong CP problem. 

To explain the motivation of this paper, we begin by briefly discussing Coleman's 
solution to the cosmological constant problem (see Section |2] for the details of the 
derivation of the following equations). We start with the path integral of the Euclidean 
gravity. If we take microscopic wormhole configurations into account, the following 
interaction AS is induced in addition to the original action; 

AS = 5^(a, + al)a I d'x^O,, (1) 

where the annihilation and creation operators of the type i babyuniverse. 

Then, the partition function of the parent universe is given by an integral over the 
eigenvalues of + aj. 

For example, if we focus on the identity operator (9 = 1, the partition function 
becomes 

^universe = y 'DgdAexp{- J d'^Xy/g{R + 2A)), 

where the wormhole effect results in the integration over A. The path integral over 
the metric g can be approximated by a 4-sphere solution, whose action is proportional 
to Therefore we have 

^universe ^ J d A 6^ , 

and the integrand has a strong peak at A ~ 0. Furthermore, if we consider the multi- 
verse, in which universes are connected each other through baby universes (see Fig{T]), 



^In [3], Banks also discussed the effect of bi- local interaction. In this paper, we mainly follow 
Coleman's argument. 

^Although [T2HI4] also studied the wormhole effect in the Lorentzian gravity, our mechanism is 
different from the previous work as we will discuss in Section [5731 
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the above integrand is replaced to exp(exp(-i-)), and the peak gets stronger. Based 
on this argument, Coleman claimed that the cosmological constant problem could be 
solved by the wormholes. 




Figure 1: A sketch of an example of the Euclidean multiverse. Parent universes are 
interacting through baby universes. 



What does this argument imply to Lorentzian spacetime? Naively, the 4-sphere 
solution is interpreted as a bounce solution. Therefore, the exponential of the action, 
ei, is expected to give the amplitude of a universe tunneling form nothing to the size 
of the 4-sphere (see Figj2]). However, if we computes the tunneling amplitude directly 



a = 




Figure 2: The 4-sphere solution can be interpreted as a foliation of 3-spheres whose 
radius expands from zero to -4t and then shrinks to zero. 



by the WKB method as Vilenkin did [15], we obtain a factor e~^, instead of ca. In 
this sense, the physical meaning of the 4-sphere solution is not clear, and whether or 
not Coleman's mechanism works in the physical Lorentzian spacetime is doubtful. 

In this paper, in order to clarify this point, we study the wave function of the 
Lorentzian multiverse consisting of infinitely many parent universes which are inter- 
acting with each other via wormholes |T6]. We will show that the density matrix 
of one universe has a strong peak in the space of the coupling constants induced by 
the wormholes. This indicates that "the big fix" indeed occurs, that is, the coupling 
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constants are determined dynamically by the quantum gravity. In particular, the mul- 
tiverse wave function predicts that the cosmological constant in the far future becomes 
extremely small. We will also find that the wormhole effect fixes the other coupling 
constants such as the Higgs parameters and the strong CP phase in the standard 
model. 

This paper is organized as follows. In Section HI we review the derivation of the 
effective action ([T]) and obtain its Lorentzian counterpart via a Wick rotation (see FigJH 
which is the Lorentzian version of FigJT]). We see that because of the wormholes the 
wave function of the parent universes becomes a superposition of states with various 
values of the coupling constants {Xi}. 

In Section [3l for the fixed coupling constants |Aj}, we calculate the wave function 
of a parent universe (j)E=o{z), where z is the sizqj of the universe, by using the WKB 
approximation. We assume that the parent universes have the topology of 5''^, and use 
the superminispace approximation for each of them 

ds^ = a\-N{tfde + a{tfdnl), (2) 

where dVL\ is the metric of unit S^. We also assume that they are created from nothing 
at a small size e via some tunneling process. Then the wave function of each parent 
universe is given by 

1 r / / 

(t)E=o{z) = sin( / kE=o{z')dz + a), (3) 

y/ii/2^Jz^kE=Q[z) Jo 

where £" = represents the so-called Hamiltonian constraint, which we will discuss 
later, and kE=o is defined by 

fc|=o(^) = -2f/(z) (4) 
= 9A-^K+ ^ + _ 5 

Here A, Mmatt, Srad are the cosmological constant, the amounts of matter and radiation, 
respectively. In principle, they are determined by examining the time evolution of the 
universe, once its initial condition at 2; = e and the coupling constants {Xi} are given. 
In this sense, they depend on the coupling constants {Xi} as well as on time, or z. 

A = A({Ai}, Z), Mmatt = Mmatt{{Xi}, z) , Srad = Srad{{Xi}, z). (6) 



^Strictly speaking, z = a"^/9 has a dimension of volume. However, for the sake of simpHcity, we 
call it "size" . 
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For example, if some matter decays into radiation at some Srad increases at this 
point. The factor . in (I3]) behaves hke for large z and plays an important 

role for our mechanism. 

In Section HJ we construct the wave function of the multiverse. The N-universe 
state |$Ar) is obtained by taking a tensor product of N universes and superposing over 

$iv(zi, ■ ■ ■ ,2;iv) ~ / rfX yu^0ij=o(-2i) ■ ■ ■0i?=o(2;iv) (S) w(\)\\), (7) 



— * — * 

where A represents the set of induced coupling constants {Aj}. |A) is the eigenstate of 
ai + a\ with eigenvalue {Aj}, and w(A) is the initial wave function of the baby universes, 
/i is the probability amplitude of creating one universe. Then, the multiverse state can 
be written as 

oo 

|0mn/ii) = ^ |$iv), (8) 

where \^n) is the N-universe state whose ^-representation is given by ([7]). Then the 
density matrix of our universe is obtained by tracing out the other universes. The 
summation over N results in an exponential, and we have 

p{z',z)(x j dX \w{X)\^\p\^ (pE=o{z'y(pE=o{z) X exp(^j dz''\fi(j)E=o{z'')\^^. (9) 

In Section |5l we try to fix the cosmological constant A by examining the A depen- 
dence of the above density matrix. If A < 0, 4>e=o is exponentially suppressed for large 
z, and the exponent in the RHS of the density matrix (Q takes some finite value. On 
the other hand, if A > 0, it is calculated as 

because the leading behavior of the momentum for large z is given by k'^^Q = 9A + • ■ ■ . 
Since this integral is logarithmically divergent, we introduce an infrared cutoff ziji for 
z, so that the above integral becomes 

(11) 

\/9A 



Thus we find that the integrand of (Q has an infinitely strong peak at A = 0, which 
means that the cosmological constant in the far future is automatically tuned to zero. 
Although we can not specify the origin of zjr at this stage, it is natural to consider 
that a sort of infrared cutoff should appear in any constructive definitions of quantum 
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gravity. For example, in the dynamical triangulation [T7], the number of simplexes 
corresponds to the infrared cutoff, or in matrix models, in which space-times emerge 
dynamically, it is provided by the size of the matrices. 

However, there is a subtlety here. There is a critical value of A = such that for 
A < Acr, a classically forbidden region, fc|;=o < 0, appears in z-space (see Figure [7j), 
and a tunneling suppression factor should be multiplied to ffTT]) . Thus, for fixed Srad 
and Mmatt, the density matrix becomes maximum when A=Acr- For example, if we 
assume the radiation dominated universe and set Mmatt = 0, we have Acr = ^/Srad, 
and the cosmological constant is fixed at 

A=A,, = 1/Srad. (12) 

Once it is done, ffTUl) becomes 



l/^l 



|2 



/ dz" OC Srad log^/R. (13) 

J ZKe=0 



Recalling that Srad also depends on the induced coupling constants the above 

equation shows that {Aj} are fixed at the values where Srad becomes maximum. There- 
fore, the value of A is given by 

A~l/max5,,rf(A). (14) 

A 

Since Srad is proportional to the volume of the universe, if the universe is sufficiently 
large, Srad is large and A is close to zero. 

To summarize, the wormhole effect makes the wave function of the multiverse a 
superposition of various values of coupling constants, but they are fixed in such a way 
that the radiation in the far future is maximized. We call it the big fix following 
Coleman. In particular, the cosmological constant is fixed as its value in the far future 
becomes almost vanishing. 

We can give an intuitive interpretation of the above mechanism. The exponent in 
the density matrix (Q turns out to be the time that it takes for the universe to expand 
from the size e to zjji. To see this we rewrite it as 

f dz \(j)E=o{z)\^ = [ dz ■ \ , = [ dt, (15) 
J Je zkE=o{z) J 

where we have used the classical relation k ~ z/z. Thus, the exponent is nothing 
but the lifetime of the universe. Naively, smaller A is favored because then the 
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universe expands slowly (see Figure [9]). However, for A < Acr, the universe bounces 
back to a small size in a finite time. Therefore, the lifetime of the universe becomes 
maximum when A = Acr- We note that the enhancement arises from the large z 
region z ~ zj^, where the universe can be described by the classical mechanics, which 
justifies treating the matter and radiation classically as in ([5]). On the other hand, 
the quantum mechanical nature of the wormholes reflects in superposing the states 
with various {Aj}. In Section [531 we compare our mechanism of the big fix with the 
previous works by other authors. 

In Section [HI as an illustration of the big fix, we consider the parameters in the 
Higgs sector in the standard model, that is, the VEV Vh and the quartic coupling 
constant X^. We assume that the other coupling constants are fixed to their observed 
values. We consider the case that Srad in the far future consists of the decay products 
of protons. Then, we can show that Srad is maximized when N^rripTp is maximized (see 
around f[M|) ). where N^, mp and Tp are the total baryon number before the decay, the 
proton mass and the proton lifetime, respectively. Naively, this seems to be maximized 
when rrip = mp{vfi) is minimized because in the usual GUT we have 

Tp oc rUp^. (16) 

Then, the wormhole mechanism seems to select out Vh = because the proton mass 
rUp depends on monotonically as follows 

rupivh) = M^°) + 3 X mu,d{vh), (17) 

where Mp^^ is the proton mass in the chiral limit, and m„ ,i is the current quark mass, 
which is proportional to Vh- However, the mass of the decay products also depends 
on Vh, and as we will show, it is in fact possible that rripTp becomes minimum at some 
nonzero value of rjiu^divh)- 

Assuming that the Higgs VEV is fixed at the observed value, i.e. 246 GeV, we next 
consider the Higgs mass. A/i-dependence enters into the above combination N^nipTp 
through the sphaleron process if we assume the leptogenesis. Smaller makes the 
sphaleron process happen more frequently and produces more baryons A*";,. Combining 
this with the fact that the stability of the potential requires a lower bound on Xh, we 
can deduce that the smallest possible value of Xh is chosen by the big fix. This means 
that the physical Higgs mass should be at its lower bound, that is, around 140 ± 20 
GeV [18]. 
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We then consider the strong CP problem. We analyze how the combination N^m^Tp 
depends on 9qcd, and find that it becomes maximum at Oqcd = 0, which means Oqcd 
is fixed to zero by the big fix. 

In Section [TJ we study universes with other topologies than S^. So far, we have 
assumed that all the parent universe have the topology of 5''^. If we allow universes 
with various topologies to emerge, we must sum over them in the multiverse wave 
function. Then, the density matrix is modified to 

p{z',z) OC j C/A W(A)VP (pE=o{z')*(pE=o{z) X 6XP^^ J '^/|/ia"0s=o(^")r^'(18) 

where a labels the topology of the universe, and /i^ is the probability amplitude 
to create such universe. Thus, the exponent of the density matrix is the sum of 
contributions from various topologies. We repeat the same analysis as for the other 
topologies, and find that the fiat universes have the largest contribution. In 

this case, the vanishing of the asymptotic cosmological constant is still valid, while the 
analysis of the big fix is modified rather drastically. 



2 Effect of Baby Universes 

We first review Coleman's argument on the effect of the baby universes p](see also ^). 
We start from the Euclidean Einstein gravity with a bare cosmological constant Ag, 

j Vgexp{-SE) = j Vg^^expi- j (fx^{R + 2K^)). 

A Planck-size wormhole configuration effectively adds to the partition function the 
following bi-local interactions (see Figure J3]), 

jvg]^ Q,e-25- j d'xd'y^^./M0\x)0^{y) exp(-5^), (19) 

where the repeated indices i,j are contracted. Cjj are some constants, and 2S'„h is 
the action of the wormhole. Summing over the number of wormholes amounts to the 
factor 

exp( -e"^"^"'' 



d'xd'y^^)^MO\x)0^{y)]. 



^Since we assume the universes are spatially compact, the topology of flat universe is actually 
torus. 
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Figure 3: A wormhole induces bi- local interactions at its legs. 



By introducing auxiliary variables Aj, the bi- local interactions can be rewritten as local 
interactions as follows, 

j l[dX, exp(^-e-^-'^A, j d^x^^)0\x) ~]^\(r^X^, (20) 

where (P^ is the inverse of the matrix cij. For example, the identity operator 0^{x) = 1 
[i = 1) shifts the bare cosmological constant Aq linearly; Aq — t- Aq + e^'^^'^Aj, which 
becomes a variable to be integrated over. 

Alternatively, we can express the wormhole effect by using the following Lagrangian 

SoS = Se + e-^"-^ 5^(4 + a,) f d^x,/^0\x), (21) 

i 

where we have introduced pairs of operators ai and a\ satisfying [ai,a^] = Qj, which 
can be interpreted as the creation/annihilation operators of the baby universe of type 
i. To understand this formula, one considers an amplitude between the initial and 
final state both with no baby universe {Vt\ exp(^e~'^"'^ + '^I) / d^x^/gO^) \ fl). By 

using the Baker-Campbell-Hausdorff formula, it is easy to show that this amplitude 
recovers Eqn. fl20l) . Although fl20l) and fl2T]) are equivalent, fl2T|) is more convenient to 
construct the wave function of the universe. 

Finally, we obtain the Lorentzian counterpart by the Wick rotation, 

S = J d^x^-g{x){R - 2Ao) - e^^""^ ^(aj + a^) j d^x^ -g{x)Oi{x). (22) 
We use this action to study the naturalness problem. 

3 Wave Function of the Universe 

In this section, we forget about the wormhole effect for a while, and consider the wave 
function of a parent universe for the fixed coupling constants Aj. We quantize the 
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system of the mini-superspace via path integral, and determine the wave function by 
the WKB method. However, as we will discuss in section |5l the whole picture about 
the big fix does not depend on these approximations, but holds quite generally. 



3.1 Wave Function of a Parent Universe 

We start from the Einstein-Hilbert action, 

j Vg^^exp{iSA) = j Vg^^exp{ij^-^ j d^x^/ -g{x){R - 2A)). 

We will consider the homogeneous, isotropic and spatially compact universe: 

ds' = a\-N{trdt' + a{trds%^,,J, (23) 

where = and dsl^^^^^^i is the metric on the spatial hypersurface, which has a 
constant curvature Ka = 1, 0, —1, depending on its topology 
Substituting the metric (123|) . the action becomes 

Sa = -\Jdt N[adyN^ - {K^a - Aa^)] , 

where we have written by the same symbol "A", which is the dimensionless 
cosmological constant. In terms of z{t) := ^^"5 it can be expressed as 

= j dt N[zyzN^ - (ir„(9^)^/3 - 9A^)]. 

The momentum conjugate to z is given by Pz = —z/zN, and the Lagrangian can 
be rewritten in the canonical form. 

La = PzZ - NUa, 

where 

1 9^/3 q 

HAip^.z) := z[--pI-U{z)1 where f/(^) := ^ir, - -a. (24) 

To describe a more realistic universe, we need to consider the energy densities of 
various fields. Then, instead of ( IMl) . the potential U is replaced to 

'Ha{Pz,z) := z[-^pI-U{z)], 

q1/3 9 A/f 9Q 

With 2Uiz) = -9A + ^K^-^^^-^, (25) 



^The spatial topology of the universe is torus and sphere for =0,-1 respectively. However, 
there are many topologies for = — 1- 
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where the last two terms represent the radiation and matter energy , respectively, and 
the associated powers of z are determined by their scaling behavior, pmatt. oc 
and prad. oc a{i)~^. We note that the coefficients depend on z and A^; 

A = A({A,}; Z), Mmatt = M^att{{\^}■, Z), Srad = Srad{{Xi}; Z). (26) 

In principle, they can be determined by solving the time evolution of the theory with 
coupling constants {Aj}, if the initial condition of the universe is completely specified. 
For example, A changes its value at the end of the inflation, and a portion of Mmatt 
may convert to S^-ad when some matter decays into radiation. 

To quantize this system via path integral, we take the following metric on the 
configuration space 



\\5g,.\\' = I d'x^gg,,g,^5g^P5g^^ oc / dt^UsNf + Na{5af), (27) 



which is invariant under the general coordinate transformation, and leads to the volume 
form of the functional integral 

n a^5N5a oc mN5z := [dN][dz]. (28) 

Collecting these results, we find that the universe is described by the following path 
integral, 

[dN][dz][dp^]exp{i / dt{p,z - NHa)), (29) 



where is given by 

In the rest of this section, we will determine the wave function of the universe, 
assuming that it initially has a small size e (see Fi^, The amplitude between z = e 




z = e 



Figure 4: The path integral ( l3T|l is defined as a sum over all histories connecting two 
geometries. 
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and z = z is given by the following path integra 

{z\e-'^\e)= [ [dp,][dz][dN]exp{i [ dt {p,z - N{t)V,K)). (30) 
J Jt=o 

z{0)=e, z(l)=z 

By choosing the gauge such that N(t) is a constant T, the path integral of N(t) is 
reduced to the ordinary integral over — oo < T < oqj, 

dT / [dp^] [dz] expii / dt {p,z - TUa) I 

J-oo J \ Jt=0 / 

z{0)=e, z{l)=z 

/oo 
rfT(z|e-*™^| e) 
-oo 

= Cx {z\6{'HA)\e) 

= Cx{z\5{HK)(^j rfE|0s)(0s|^|e). (31) 

From the first line to the second line, viewing T'Ha as the Hamiltonian, we have used 
the ordinary relation between the operator formalism and the path integral one, and 
C is some constant. In the final line, we have inserted the complete set {I^e) } defined 
by 

{ct>'E\<pE) = 6{E-E'), (32a) 
UA\ct>E) = E\cl>E). (32b) 

Therefore, by using (j)E{z) = {z\(j)E)i the amplitude can be expressed as 

C X 0^,=o(e)feo(^). (33) 

In other words, the quantum state of the universe that emerged with size e is given by 

Cxcl)%^^{t)\(l>E=o). (34) 

We can calculate 0_b(-2) in the canonical quantization formalism. By replacing 
Pz — 7- —id/dz in the Hamiltonian (1251) . Eqn.( l32bl) becomes 

^^{\^2 - U{z))^z 4>e{z) = EM^)- (35) 



®This analysis is similar to that of [14] . 

^To be precise, we should integrate only positive T if we fix the time-ordering of the surface St=o 
and St=i as in Fig|4l However, we take the integration range as — oo < T < oo to obtain the 
well-known Wheeler-Dewitt equation in the path integral formalism. This procedure corresponds to 
summing over the ordering of the two surfaces too. 
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Note that for E = this leads to the Wheeler-DeWitt equation. However, we need to 
solve this equation for general E since we should determine the normalization constant 
of the wavefunction according to fl32ap . We rewrite ( l35l) as 



dz"" 



- kl{z))yGMz) = 0, (36) 



where 



kUz) ^ -2Uiz)-^ 



Q._9^^ 2MmaU , 2Srad 2E 
^2/3^"+ ^ + ^4/3 ^ ' 

and apply the WKB method to the function ^/z(f)E{z)■ The solution in the classically 
allowed region, k'^^Q^z) > 0, is given by a linear combination of 

1 r I I 

(pE=oiz)= ^ ^ r- ==exp(±i / dz'kE=o{z')), (37) 



where the normalization is determined by f l32ap (see Appendix A ). 

We need to specify the boundary condition to determine the solution completely. 
As a simple example, if we require 0£;(O) = Ojlj we have 

0i;=o(^) = ^— ^ , , ^ sin( / dz'kE=o{z')). (38) 



V^^Vz^/k^I^ 

However, we do not need the details of the solution in the following sections. 

4 Multiverse Wavefunction and Density Matrix of 
our Universe 

In this section, we construct the multiverse wave function assuming that all the parent 
universes have the topology of 5^. Here, we mean by the word "multiverse" the state 
with an indefinite number of universes. We then calculate the density matrix of one 
universe, which is essentially what we observe in our universe. 



^°The boundary condition would be more complicated because the behavior in z < e is determined 
by the dynamics near singularity. 
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Figure 5: A sketch of an example of the multiverse. Each parent universe emerges 
with a small size e by a tunneling process. In this example, the initial state has no 
baby universes and the final state has two baby universes. 



4.1 Wave Function of the Multiverse 

Usually, the universes which are not connected with ours are irrelevant for us, since 
they have no effect on our observation. However, when we take the wormholes into 
account, these universes interact through them, and all the universes become to have 
the same coupling constants which should be integrated in the path integral. 

In order to construct the quantum state of the multiverse, we need to specify the 
initial state of the baby universes, which can be expressed as a superposition of the 
eigenstates of the operators ai + a], 

{ai + al)\\) = \,\\), (39) 

where we have denoted the set of coupling constants {Aj} by A. For example, if there 
are initially no baby universes as in Fig. \5\ the state is given by J Yli dXi e~^^^''''^^^^\X) := 
where \Q) satisfies ai\Q) = OO In general, there may be many baby universes 
initially (see Figurej6]), and the state can be written as jYlid-K w{X)\X), where w is 
a function of A. 

To write down the multiverse state, we also need the probability amplitude of a 
universe emerging from nothing, which we denote by fio in analogy of the chemical 
potential. Here we assume that all universes are created at the size e. Together with 
the factor in (IMl) , the weight of each universe is given by, 

/i:=/ioxCx0*^^o(e). (40) 



"'^"'^It might be helpful to regard a + as the position operator \/2x of a harmonic oscillator, and 
recall the ground state of the system |0) can be written as / dx e~^^^^ in the cc-representation. 
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Figure 6: A sketch of an example of the multiverse. In this case, the initial state has 
some baby universes. 



A crucial fact is that /i does not depend on A strongly. This is because 0^^Q(e) is a 
smooth function of A as is seen from f p8|) . and C arising from the path measure should 
have nothing to do with Aj. 

Then, the multiverse wave function can be written as 



\<Praulu) = 1*^^) (41) 

N=0 

where \^n) stands for the A^- universe state, whose wave function is given by 

^Nizi,- ■ ■ ,zn) = dX X (f)E=o{zi)(f)E=oiz2) ■ ■ ■(l)E=oizN) w{X)\X), (42) 



where 

dX = Y[dXi. (43) 

i 

4.2 Density Matrix of Our Universe 

We now can obtain the density matrix of our universe by tracing out the other universes 
and the baby universes, namely A. Using ( |42|) . we can calculate it as 



f dz^ 

p{z',z) = Y <i>%^^{z',zi,--- ,zn)'^n+i{z,zi,--- ,Zn) 

= dXwiXf\fi\^(l)E=o{z'y(i)E=oiz)x I / dz"\fx^E=o{z")\ 



N=0 



dX w(A)Vr (l)E=o{z'y(f)E=o{z) xexpl^J dz \i2(j)E=o{z )pj, (44) 

where z and z' are the size of our universe. We note that the above integrand depends 
on {Xi} through the wave function 4>e=o- 
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5 Vanishing Cosmological Constant 



In this section, we show that the integrand in ( l44l) has a strong peak at a point in the 
{Aj} space where the cosmological constant A = A({Aj}) becomes very small, which 
means the cosmological constant problem is automatically solved. We also discuss the 
possibility of the big fix. 

5.1 Evaluation of the Density Matrix 

In this subsection, we examine how the exponent in the density matrix (jH]), 

dz"\n(f)E=o{z")\'^, (45) 





depends on A. 

First we sketch the potential U{z) in (125|) . Again we assume that all the universes 
have the topology of S*^ {K = 1), so that U{z) is given by 

q1/3 9 A/f 9Q 

For large z, the leading term is the cosmological constant A, and the next leading term 
is the curvature term. We note that only the curvature term is positive, and U{z) has 
a maximum at one point z = z^,, 

U\z,) = 0. (47) 
As we vary A with Mmatt and Srad kept fixed, U (z) changes as in Fig [71 There is 



a critical value Acr at which the maximum becomes zero (see Fig 7(d)); 



U{z,)\a=a^^=0. (48) 

Note that if A = Acr, three contributions to U{z), the cosmological constant, curvature 
and energy density coming from matter and radiation, are comparable around z z^,. 
The precise values of z^, and A^r depend on the history of the universe. If all the matter 
decay into radiation hj z = z^:, we have Mmatt = 0, and Acr is given by 

g^3/2 gi/3 

z^ = — —, 9 Acr = — • (for radiation dominated) (49) 

3 8 Srad 

On the other hand, if the matter dominates around z^, they are given by 



•iMmatt QA _ l 

~9V^' ~ 3M? 



9^cr = 77772 — • (^°^ matter dominated) (50) 

matt 
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(a) A < (b) A = (c) < A < A^r (d) A = A^^ (e) A > A, 



Figure 7: As we vary A from zero to Kcr ~ -5 — , the region where the wave function 
takes the tunnehng suppression becomes shorter. For A > A^r, there is no suppression. 



Now we can examine the behavior of (I)e=q{z) in the large-z region, and evaluate 
the integral (H5|) . If A < 0, the wave function damps exponentially, and (H5|) is finite 



(see Fig 7(a)). On the other hand, if A > 0, the wave function does not damp for 
sufficiently large z, and f HS]) is divergent. Thus, if we introduce a cutoff for large z, as 
we will do below, (H5|) takes the maximum for some positive value of A. 

Furthermore, if A > Acr, all the region of z is classically allowed, and we can 
reliably use the WKB solution 

'^''='^'^ ^ /zh (zV ^^^^ 

\/ ZKe=0\Z) 



which becomes larger as the momentum fc_B=o = a/— 2C7 becomes smaller. Thus, for 
A > Acj., the wave function becomes the largest when A = K^.,.- On the other hand, 
if < A < Acr, there is a forbidden region, which suppresses the wave function. The 
suppression is stronger for smaller A because the forbidden region becomes larger as 
we decrease A. Thus, we find that fHSj) takes its maximum value at 

A = K^r- (52) 

Next we discuss how the maximum value of ( H5|) is determined by the amount of 
radiation Srad or matter Mmatt- If we set A = Acr, using fl5T|) we have 

dz \fMj)E=o{z )P~ l/Up / dz — — — . (53) 



Z^/Ar 



'0 ^0 

Since this is divergent, we introduce an infrared cutoff zm and replace z = oo with 
z = ziR. Then the above integral becomes 

1 1 

dz ~ -j= log ziR, (54) 

Z\ YiVcj. 
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The cutoff zjR should be explained from a microscopic theory of gravity such as string 
theory. For example, in the IIB matrix model space-times emerge dynamically from 
the matrix degrees of freedom, and an infrared cutoff appears effectively, which is 
proportional to some power of the matrix size [T9]42T| . 

If we consider the case of fH9l) . where the curvature term balances with the radi- 



ation, flM|) is proportional to ySVod log zm, and the integrand of the density matrix 
(HH) behaves as 

^ ^ (55) 



exp|^COrast. X a/ Srad log Zl]^ , 



which has an infinitely strong peak at a point in the {Aj} space where Srad becomes 
maximum. Here, we have assumed that does not have a strong dependence on 
{Aj} because it is determined by the microscopic dynamics of smaller scales than the 
wormholes. Thus we have seen that all the couplings {Aj} are fixed in such a way 
that Srad is maximized. We call it the hig fix following Coleman. In the original 
Coleman's argument the enhancement comes from the action itself, or equivalently, 
the exponential factor in the wave function fl37|) . while it comes from the prefactor in 
our case. We will discuss this meaning in the next subsection. We also note that the 
big fix applies only to the couplings that are induced by the wormholes. In particular, 
the cosmological constant is given by 

A = 1/max ^^.^^(A), (56) 

A 

which is very closed to zeroo We note that A and Srad appearing above should be 
regarded as their values ai z = z^,. 

In the other case (!50l) . where the curvature term balances with the matter, we have 
Acr ~ ^mait) instead of fl55l) we have 

exp (^onst. X Mmatt log zi^ . (57) 

This time, the coupling constants {Aj} are fixed such that Mmatt sX z = z^, is maxi- 
mized, and the cosmological constant is given by 

A = l/maxML,,(A). (58) 

A 



^"^Srad means the amount of the radiation in the whole S''^-universe, rather than that in the portion 
we can observe. Thus, if the whole universe is large enough, Srad is extremely large. 
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In the above mechanism, the curvature term becomes comparable to the cosmo- 
logical constant around z = z^.. On the other hand, observational cosmology tells that 
the former is much smaller than the latter already in the present universe. Therefore, 
in order for the scenario to work, the cosmological constant needs to decrease as a 
function of time by some mechanism such as quintessence models. Then the above 
argument claims that its asymptotic value is very small. 

5.2 Interpretation of Enhancement at A = 

In this subsection, we provide an intuitive explanation of the enhancement at A = k^r 
in (144 p . We also argue that our mechanism works beyond the minisuperspace and the 
WKB approximation. 

First, we recall that the enhancement of the density matrix comes from the expo- 
nent in (jS]), 

y dz |0s=o(^)r, (59) 
and we have evaluated it by using the WKB solution 

0s=o(^) ~ , / , (60) 
Classically fc_E=o(2;) is the conjugate momentum of z, 

kE=oiz) ~ z/z. (61) 

Thus, (159|1 can be written as 

r r^iR 1 r^iR 

/ dz |0^=o(^)P = / dz = / -, (62) 

J Je zkE=o[z) J, Z 

which is nothing but the time it takes for the universe to grow from the size e to 
zjR. Since we have imposed the cutoff Zir on the size of the universe, a universe with 
the size larger than zir does not existo. Thus, ( 162|) can be interpreted as the time 
duration in which the universe exists. We call it the lifetime of the universe, for 
simplicity. 



"'^'^ Although we have not specified the infrared cutoff precisely, we can simply imagine that when a 
universe reaches the size zm, it ceases to exist , or it bounces back and starts shrinking towards the 
size e. 
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In fact, we can verify this interpretation without relying on the WKB approxima- 
tion. We recall the normalization of the wave function 



b'E\<l>E) = 6{E - E') 



(63) 



which leads to 

J dz |0s=o(^)r~5(O). (64) 

As is usually done in the derivation of Fermi's golden rule, (5(0) is regarded as the 
total interval of time, which in our case is naturally interpreted as the duration of the 
universe. 

Therefore, what the big fix does is to make the lifetime of the universe as long 
as possible. Based on this interpretation, we can reproduce the results obtained in 
the last subsection. First we note that, for A < Acr, the universe cannot reach to zm 



because of the potential barrier (see Fig, 8(a) ), and collapses back to the size e and then 



disappears in finite time (see Fig|9](a))o So we concentrate on the case A > Ac,-. As 
we vary A, the depth of the potential changes as in Figure [HI The shallower potential 
gives the longer lifetime, and thus the lifetime becomes maximum at A = K^r (see 
Figil^b) and (c)). 




-9A 




-9A 




(a) A < Ac 



(b) A ~ Ac 



(c) A > Ac 



Figure 8: The classical motion is shown for each A. For A > Acr, the universe expands 
to zjFt rapidly, and the lifetime is short. For A ~ Acr, it takes long time to reach zjr, 
that is, the lifetime is long. 



Before closing this subsection, we emphasize the general validity of our mechanism. 
So far, we have used the mini-superspace approximation, in which only the size of the 

^^Quantum mechanically, the universe can reach to zir after tunneling for < A < Acr, but 
because of the tunneling suppression such A does not contribute much, as we have discussed in the 
last subsection. 



19 



ZIR 








(a) A < A, 



cr 



(b) A ~ A, 



cr 



(c) A > A 



-cr 



Figure 9: For A > Acr, as A varies to Acr, the universe takes more time to expand to 
the size zir, and the "hfetime" becomes longer. For A < A^,., the universe bounces 
back to the size zero before reaching zir. 

universe is considered, and the other degrees of freedom such as various fields and 
inhomogeneous fluctuations of the metric are ignored. If we take those degrees of 
freedom into account, the quantum state of the universe is described not only by z, 
but also by the other degrees gi, and ( l59l) is replaced by 



However, if the quantum state of g^'s is normalized to 1, the integration over g^'s 
leaves the same integral as the mini-superspace, and again we have 5(0). Therefore, 
we can say quite generally that the exponent of the density matrix is the lifetime of 
the universe. Furthermore, the integral ( l65i) is controlled by large values of z, where 
the evolution of the universe is completely classical. In such late time, the effect 
of the other degrees of freedom such as gravitons, photons, and protons is simply 
represented by the energy density in the potential (H^ . which justifies the analysis we 
have employed above. 

5.3 Comparison with Euclidean and Other Lorentzian Ap- 
proaches 

In this subsection we discuss the difficulty of the Euclidean gravity, and explain how 
our mechanism is different from the original Coleman's or the subsequent Lorentzian 
approaches. 




(65) 
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5.3.1 Wrong Sign Hamiltonian 

In order to clarify the problem, we start with a Hamiltonian 

2 

H^ = -j-Viq), (66) 

which is the minus of the normal Hamiltonian 

H^ = ^ + V{q), (67) 

where p is the canonical momentum of q and V{q) is a potential. Since the Schrodinger 
equation 

i^^{q,t) = H^{q,t) (68) 

for ( |66l) and (1671) are simply related by the complex conjugate, they should describe 
the same physics. In particular, the tunneling phenomena are the same: When we 
consider a tunneling process, the wave function should decrease in the direction of the 
penetration, and the tunneling is exponentially suppressed for both cases. 

Next we discuss the Wick rotation of the wrong sign Hamiltonian (|66|) . Usually, 
for the right sign Hamiltonian f l67|) . we rotate the time t = —ite so that the 

transition amplitude 

(g'|e-*^+*|g) = (g'|e-^+"^ |g) = (g'|e-"^('^+^(''))|g) (69) 

is well defined. Note that the rotation in the opposite direction t = ite does not work 
because of the bad large-momentum behavior. On the other hand, for the wrong sign 
Hamiltonian H_, we should take t = ite 

(g'|e-'^-*|g) = (g'|e^-"^|g) = (g^-^^^'^+^^'^^^lg), (70) 

and t = —ite does not work. 

Obviously, (169|) and (!70|) are the same, and thus the equivalence of the two systems 
can be seen also in the Euclidean framework. However, the Wick rotation should be 
done in such a direction that the transition amplitude is well defined. In other words, 
if one applied the naive Wick rotation t = —ite to the wrong sign Hamiltonian H_, 
one would have physically unreasonable results. 

As an example, we consider the Hamiltonian H_ with V{q) = \{q^ — q^Y and the 
transition amplitude 

(g' = +go|e-*^-*|g = -go). (71) 
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If we perform the correct Wick rotation t = ir^, the amphtude is given by the ordinary 
Euchdean path integral as is seen from fl70|) : 

(g' = +go|e^-^-|g = -go) = J Vqexp(^- j dr (^^(drq)' + V^^ (72) 

The one-instanton solution connecting q = —go to g = +qo contributes as 

{q' = +go|e-^-"^|g = -go) ~ Ce-^^[«='l + ■ ■ ■ , (73) 

where SE[qci] is given by Se = J^g^ dq^y2V{q) = ^^q^y/X. This is consistent with 
the suppression of the tunneling. On the other hand, if we perform the wrong Wick 
rotation t = —ite, the amplitude is formally given by an Euclidean path integral for 
unbounded action 

(g' = +go|e-^-^-|g = -go) = J Vqexp(^j dr (^^{drq)' + V^^ (74) 

Although this path integral is ill-defined, if we naively evaluate it by using the instanton 
solution qd, we have a wrong answer 

(q' = +q^\e-^-^E |g = _ Ce'^^^'^^'^ + ■■■ . (75) 

This would indicate that the tunneling is not suppressed but enhanced exponentially. 
However, as we have discussed above, we do not regard it as true. 

5.3.2 Case of the Quantum Gravity 

We now turn to the case of quantum gravity, whose Hamiltonian is schematically given 
by 

^ = ^[-n^ + /(a)nLns] + --- , (76) 

where Utrans stands for the canonical momentum of transverse modes of the metric, 
and /(a) is a positive function of a. We note that the signs in front of Ua and Utrans 
are opposite. The dots represent various matter and gauge fields, which have the 
same sign as the transverse modes. Thus, if we perform the standard Wick rotation 
t = —IT to make the transverse and matter sectors well-defined, we lose control of the 
fluctuation of the conformal mode. On the other hand, if we take t = +iT to avoid it, 
then the transverse and matter sectors become divergent. Thus, the time axis cannot 
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be rotated in any direction, and the Euclidean gravity obtained by a simple Wick 
rotation is problematic]^ 

In order to clarify the origin of the confusions about the Euclidean gravity, we 
consider the tunneling nucleation of the initial universe. The Hamiltonian in the mini- 
superspace is given by 

Hgrav = -a^ + Pvaca^), (77) 

where p^ac is the vacuum energy of the universe in the planck epoch. Classically, 
the evolution of a{t) is given by solving Hgrav = 0, and in quantum mechanics, it is 
promoted to the constraint on the wavefunction of the universe, 

(^-a' + p.aca^)^(a) = 0. (78) 

As Vilenkin showed by using the WKB analysis [15j, the tunneling probability V from 
a = to a = 1/ y/pvac is given by 

2 

VwKB oc e ^pvac , (79) 

This result can be obtained in the Euclidean formalism, if we apply the Wick rotation 
correctly, t = ite, as we have discussed in the previous subsection. Then, the bounce 
solution a(r£;) is given by 

a{TE) = ^^cos{y/p^TE), (80) 

Y Pvac 

and, for this solution, the Wick rotated action is evaluated as 



SE[a]= I dTE-a{l + {- ) - pvacCt =- , (81) 

Z \ OTe J iPv 



-'vac 



from which we obtain the tunneling probability V as 

V cx exp(-^E) = exp(--^). (82) 

'^Pvac 

We can thus recover fl7^ . and there is no enhancement as p^ac +0. 



-•^^There is some argument that the analytic continuation of the conformal mode might cure the 
problem |22] • Here we do not consider this possibility since the physical meaning of the complexified 
scale factor is not clear. 
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On the other hand, if we performed the Wick rotation in the wrong direction 
t = —ite, which is the case of the ordinary Euchdean gravity, we would obtain an 
enhancement instead of the suppression, 

P = exp(5s)=exp(-^), (83) 

which states that the bigger universe is more hkely produced via the tunnehng. It 
seems that this picture is accepted in the original Coleman's and some of the sub- 
sequent works, and used to discuss the possibility of the double exponential form 
exp(exp(^-^)) in the multiverse. However, as we have discussed, we do not accept 
this picture, and we simply trust the results of the Lorentzian gravity, in which the 
tunneling is suppressed. Therefore, we do not claim the double exponential form, and 
instead we have shown a different origin of the enhancement, which leads to (l55l) or 

m- 

5.3.3 Enhancement in the Lorentzian gravity 

Here we discuss how our enhancement mechanism is related to the probabilistic inter- 
pretation of the Wheeler-DeWitt (WDW) wave function, and compare our mechanism 
with the other authors'. 

First, we emphasize that our enhancement mechanism has a completely different 
origin from Coleman's original idea; he obtained the enhancement at A = from the 
path integral itself, which is evaluated by the 4-sphere solution as 

j Vg e"^^ ~ ei. (84) 

We think this is fake as we have discussed in the previous subsection. On the other 
hand, our enhancement mechanism has nothing to do with the value of the path 
integral. In fact, by using ( 13T|) . the amplitude of a universe emerging with Zi = e and 
terminating with = e is evaluated as 

J f/T(e|e-^^^|e) ~ VE=o(e)0*s=o(e), (85) 

which is not particularly enhanced. 

Even though the path integral itself does not have enhancement, it arises from 
the probability measure of the WDW wavefunction. In this paper, we have simply 
assumed that the absolute value squared of the wave function gives the probability 
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density |23 ||2^^1 . More specifically, the multiverse state ( HTj) is the superposition of Ad- 
verse states each of which consists of universes with sizes zi, - ■ ■ , zn-, and coupling 
constants {A}, and we interpret 

N 

\^{zi, ■ ■ ■ , ZN)\'^dXdzi ■ ■ ■ dz^ = l/ip^ Y\ \4'E=o{zi)\'^dzidX (86) 

i=l 

as the probability of findin g N universes with the sizes Zi ^ Zi + dzi{i = 1, ■ ■ ■ , A^) and 
the coupling constants {A}|^ 

Although this probability measure is a straightforward generalization of the ordi- 
nary quantum mechanics, 

(87) 

there is some criticism. If we evaluate the normalization integral 

dz\(l)E=o{z)\', (88) 



we find a divergence for large z. It essentially comes from the integral over time 
T in the path integral ( pTl) . which makes the universe a superposition of z. Thus, 
the measure fl86|) appears to correspond to the following probability measure in the 
ordinary quantum mechanics 

\(j){x,t)\^dxdt, (89) 

whose integral is obviously divergent since J \(f){x,t)\'^dx is constant in time. 

However, we adopt the probability measure ( l86l) as the most natural one. The 
reason is the following: Suppose we perform a numerical simulation of some micro- 
scopic model that realizes the emergence of the multiverse. Then, every time we make 
an observation, we find an A^-verse which consists of A^ universes with various sizes. 
Therefore, we are naturally lead to consider the ensemble of A^ universes with the 
probability fl86p . The divergence of (155]) practically does not cause any problems in 
the process of the simulation. As we have mentioned, the infrared cutoff zjji is natu- 
rally introduced, for example, as the size of the matrix when we design the spacetime 
geometry by matrixes, or the number of the simplexes in the dynamical triangulation. 

In order to understand how the enhancement arises from the measure fl5B]) . we first 
consider the single universe state. The WDW wave function of the universe ^4>e=o{z) 



^^For a review of the various interpretations, see for example |25j . 

^^Here, we omit the weight of the couphng constants, w{X) in since it does not play any 
important role in the argument. 
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represents the superposition of various universes with size z. As we have seen around 
the measure |/i</)£;=o(-2)P'^^ can be interpreted as the probabihty distribution of 
the time T that has passed after the universe emerged, 

\[i(\)E=Q{z)\^dz ^\ii\^dT. (90) 

If we integrate it over we find that each universe has the weight 

y dz\ii<^E=^{z)\^ --\[A^T^- (91) 

Here, T\ is the hfetime of the universe, which depends on the couphng constants {A}. 
Thus, our probabihty measure counts the universe with the weight \[i^T\. Similarly, 
the iV- verse state is the superposition of the states each of which consists of universes 
which were created at various times. Therefore, fl86|) is equal to 

\yS^^dT^---dT^, (92) 

where we consider that the i-th universe was created time Tj before the observation. 
If we integrate fl5^ over the sizes, the iV-verse state is counted with the weight 

^l/i|'^A, (93) 

where A^! is the symmetry factor. When we evaluate the density matrix (jS]), the 
lifetime T\ becomes exponentiated to 

exp (l/ipTA) m 

after summing over the number of the universes. Thus, our enhancement mechanism 
essentially comes from the probability measure, which counts each universe with the 
weight of the lifetime. 

We expect the big fix occurs in such a way that the lifetime is maximized. This 
point is completely different from the earlier works based on Lorentzian gravity (TTHH] , 
In particular, our mechanism has nothing to do with the initial tunneling amplitude /i. 
As we have seen from (1521) . in general depends on the various coupling constants {A} 
at the planck scale. However, what determines the lifetime of the universe is not the 
microscopic parameters themselves but the parameters at the low energy scale, such 
as the renormalized cosmological constant and the Higgs mass, and so on, and there 
is no reason that [i has a strong dependence on such low energy quantities. Thus, the 
tunneling amplitude [i does not play an important role in the big fix. 
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6 The Big Fix and the Gauge Hierarchy Problem 



One of the notorious problems of the standard model is the gauge hierarchy problem, 
which arises from the quadratic divergence of the Higgs mass. In this section, assuming 
that the wormhole effect induces the parameters of the Higgs potential, the VEV Vh 
and the quartic coupling A/i, we examine the possibility that the hierarchy problem is 
solved by the big fix. Here we take, as the low energy effective theory, the ordinary 
standard model with the proton decay at the GUT scale, and fix the gauge and the 
yukawa couplings to the observed values. In order to discuss the big fix, we need to 
know the universe in the future. Here we assume that the curvature term balances 



with the radiation after the baryons decay, which corresponds to the case of Fig 7(d) 



and Eqn. (l55l) . Such a universe is realized if, for example, the following conditions are 
satisfied: 

Condition 1. The cosmological constant is time-dependent and decreases to the 
asymptotic value before the proton decay. 

Condition 2. The lifetime of the dark matter is much shorter than that of protons. 

Condition 3. The curvature balances with the energy density while the decay 
products of baryons being relativistic. 

A comment is in order on the above conditions. If they are satisfied, the universe 
evolves like in FigJTOl Condition 1 and 2 ensure that the cosmological constant and 
the dark matter become irrelevant in the energy density, and so the baryons dominate 
the energy density. However, around the proton lifetime, the baryons decay and the 
radiation such as relativistic electrons are produced after the decay. Finally, as the 
universe expands, the leptons become non-relativistic, namely become matter, due to 



the red-shift. As we have discussed in Section 15711 we need to specify in which stage 
the curvature term becomes comparable to the energy density. Condition 3 claims 
that it happens in the third stage as is shown in FiglTUl In general, as we will discuss 
in [Appendix B[ the e-foldings of the initial inflation determines when it happens, and 



the above scenario corresponds to the values given by fll3ip . 



In Section 16.11 we discuss how the proton decay determines Srad in the far future, 
and we write it in terms of the proton mass mp, the total baryon number N^, and 
the pion mass m,r- In Section [6721 we will analyze how these quantities depend on \h 
and Vh and at what values they are fixed. In Section 16.31 we discuss the strong CP 
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Figure 10: A sketch of the potentiaL 



problem. 



6.1 Proton Decay and the Radiation 

We denote the proton decay rate by Tp, and its inverse by Tp. When the protons 
decay into radiation at some large scale factor ), the continuity of the energy 

density leads to the following relation 



AS.. 



rad 



AM, 



matt 



X a 



■pi 



(94) 

where AMmatt is the contribution of protons to Mmatt, and ASrad is the radiation 
amount produced by their decay. Because AM^att is expressed as Nf, x nip, the second 
equation of (!94l) becomes 

ASrad = mpX NhX ap. (95) 

We assume that the cosmological constant A decreases so rapidly that the universe 
is mattar-dominated in most of the time until the proton decay. Then, the Friedman 
equation (a/a)^ = '^^"f** determines the evolution of the scale factor as 



and fl95l) becomes 



rad 



,3/2 



N^mlTp. 



(96) 



(97) 



Tp can be estimated as follows. The effective interaction which induces the proton 
decay is given by enp with the coupling constant g/M^, where is the GUT scale 
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and g has the mass dimension two, g ~ ^qcd- Using the formula of the decay rate, 

^'-2^(7 ^ I §^)\M\H2^rs''H,.-P,-P.), (98) 

where Ai is the matrix element of the decay process, we have 

^ ~ ' ^^^^ 

and (|97|) becomes 

{ASra,f'ocN^g-'mJl-^) \ (100) 

6.2 The Big Fix of the Higgs Parameters 
6.2.1 The Higgs Vacuum Expectation Value Vh 

Before discussing the big fix of we note that, since we regard the yukawa couplings 
as constants in our argument, we can consider the current quark mass mu4, instead 

of Vh- 

mu,d = VhVuA- (101) 

Then, what we want to know is the value of rriu^d that maximizes the radiation amount 
ASrad in fllOOl) . Nh does not depend on Vh much if we assume the leptogenesis in which 
the baryons are mainly produced in the energy scale much higher than Vh- Therefore, 
we concentrate on the remaining quantities 

g, Trip, m^. (102) 

If » ^QCD, a simple quark counting and the dimensional analysis tell us 
that the masses and the coupling constant g are given by 

mp ~ 3 X mu,d, ~ 2 X mu,d, g oc ml^^, (103) 

which means that ASrad is a decreasing function for large m„ rf. We thus examine the 
possibility that ASrad becomes maximum at some small value of mu,d- 

We need an expression of the quantities (11021) for small rriu^d- Firstly, the proton 
and pion masses are given by 

ml = aM^^'^mu^, 

mp = M^°) + 3/3m„,rf, (104) 
29 



where Mp°^ is the proton mass in the absence of the current quark mass, and ot and /3 



are some numerical parameterslif 



Both of a and /3 are determined by the dynamics 



of massless QCD, and are independent of rriu^- Experimentally we have 

Mf ^ 910 MeV, m„,d ^5-10 MeV, (105) 
and so a takes some value around 2 < a < 4. 



On the other hand, since g has the mass dimension two, it can be expanded in 
as follows: 



g oc 



(Mf)Ml + 3/3«:^l, (106) 



where k is some parameter around < k < 

Substituting (fTOD and (fT06|) into (fTOOj) . we find that 

(AS.,.Poc(l.3&)-'(l-^^)", (108) 
where we have introduced x = -^j^- (HOSp can be expanded as 

{ASradf/^ OC 1+ (^2a-6K + 3(3^x + 0{x^) . (109) 

which indicates ASrad is a increasing function for small rriu d if 

2a-6K + 3/3>0. (110) 

If it is the case, since we have seen ASrad is decreasing for large rriu d, we can conclude 
that ASrad takes its maximum at some small itLua- 

In order to determine the concrete value of x, we need the second order term in 
(11091) . and more precise analysis of QCD is required. It would be very interesting to 
see whether or not ASrad really takes its minimum at the experimental value of x, 
g|(j < X < cjYo- If i^ works, the big fix fixes rriu^d to 5 ~ lOMeV, which implies the 
Higgs VEV to be 

Vh ~ C(lOOGeV). (Ill) 



^^Naively, proton mass is expected to be written as nip — Mp^^ + Siriu^d- However, turning on 
non-zero m„^d affects the chiral condensation. We express the total effect by the parameter /3. 
^^This range of k seems reasonable if we rewrite ()106|) as 
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6.2.2 The quartic coupling constant and the Higgs mass 

Assuming that Vh is correctly fixed at Vh ~ 246GeV, we next discuss the quartic 
coupling constant A/^, and predict the Higgs mass. 

The A /i- dependence of AS^ad is quite simple because A^ enters only A^";, in (11001) ]^° 
Since in the leptogenesis scenario most of the baryons are produced swiftly in the sym- 
metric phase, the baryon number does not depend on the Higgs parameters strongly. 
However, if we make A^ smaller, the period of the symmetric phase becomes longer. 
Thus, the number of the baryons A*";, becomes slightly increased. Therefore, A";, is a 
decreasing function of A^, and smaller A^ dominates in the density matrix (jH]). 

However, it is well known that there is a lower bound for from a stability of the 
Higgs potential. This bound corresponds to the case that the coupling \h vanishes at 



the Planck scale, or wormhole scaleO Thus, \h is fixed to this lower bound by the big 
fix. As shown in [T8j, the corresponding Higgs mass nth is around 

m?, ~ 140 ± 20 GeV. (112) 

We note that while we need some assumptions of cosmology in order to discuss w^, 
the argument of the Higgs mass seems relatively generic. (11121) can be derived only 
by assuming that the Higgs VEV is Vh — 246GeV and that the energy density of the 
universe is a decreasing function of Xh- 



6.3 Strong CP problem 

So far we have assumed that the CP violating phase 9 is vanishing since there is an 
experimental upper bound 9 < 10"^^. We can also discuss the strong CP problem by 
examining how the non-zero deviation of 9 infiuences the radiation amount ASrad in 
(HOOD. 

Fortunately, we can make an argument without knowing the specific ^-dependence 
of ASrad- The baryon number Nj, does not depend on 9 since Nf, is determined at much 
higher energy, and the remaining quantities, nip, ml, g, should respect a refiection 
symmetry due to the CP transformation: 

9 -9. (113) 

Although we have neglected w/i-dependence of Nb in the discussion of the big fix of Vh, we can 
not ignore in Nt, because Xh only appears in Nt, in (|100p . 

■^^We assume that the wormhole size is almost equal to the Planck scale. 
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Strictly speaking, the real CP transformation flips the sign of the CKM phase as well 
as 9. However, the reflection of 6 is an almost exact symmetry in the hadronic scale, 
which is much lower than the weak scale. Thus, ASrad must be an even function of 6, 
and we have only two possibilities: the point ^ = maximizes or minimizes ASrad ( 
at least locally). If the former is the case, and 6' = is the global maximum, 6 is fixed 
to zero by the big fix. It would be very interesting to examine by QCD whether it is 
really the case or not. 

We note that this argument is highly generic because it relies only on symmetry, 
and so we can still make a similar argument even when a cosmology other than that 
we assumed in this section is realized. 



7 Universes with Different Topologies 

So far we have only discussed closed universes with topology {K =1). In this 
section, we study the universe with other topologies. We first discuss the case that all 
the universes are flat {K = 0), and compute the density matrix. We find that it has 
a strong peak at A = 0. We then consider the case that all the universes are open 
{K = -l). 

Finally, we construct the density matrix in the case where various topologies are 
allowed in the multiverse state. We will find that the flat universes dominate in the 
density matrix. 

7.1 Flat Universes 

We consider the case that the multiverse consists of flat universes. For K = 0, the 
potential U (z) is given by 

2U{z) = -9A - - (114) 

If A > 0, the whole region of z is classically allowed, and the integral of the wave 
function can be evaluated as follows by using the solution ( l38ll with K = 0: 

dz\(j)E=o{z)\^ = [ dz / . . sm\ [ kE=oiz')dz')), (115) 

J 'KZkE=o[Z) J 

which is divergent and we regulate it by an infrared cutoff zm as before. Then it 
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(a) A > 



(b) A < 



Figure 11: The potential U{z) for the flat universe. The solid line is the potential, and 
the dashed line is its asymptotic value —9 A/2. The colored line represents a typical 
form of the wavefunction (1>e=o{z). 

behaves as 

/ dz\<PE=o{z)\' ^==logziR + --- . (116) 

On the other hand, if A < 0, (j)E=o{z) damps exponentially for large z, and the integral 
gives a finite value. Therefore, this region can be neglected in the density matrix. 
Then we obtain the following density matrix (jSj), 

p~/ d\ |/i|^e~^0j7;=o(-2O*0g=o(-2)exp( — x log zir), (117) 
Jo VOtt^A 

which has an infinitely strong peak at A = 0. Then, A- integration can be performed 
simply by substituting A = in the integrand, and the exponent in the density matrix 
can be written as 

dz ——— I dz ^ , =Ztti , (118) 

zk{z) J z^M^^u/z ^m^t 

where we have assumed that the universe becomes matter dominated for large z. 



7.2 Open Universes 

For K = — 1, the potential U{z) is given by 



2U{z) = -kl^,{z) = -9A - - = - (119) 



where the second term comes from the negative curvature. As in the case of the flat 
universe, this potential is always negative for A > 0, while for A < it becomes 
positive for large z. 
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If all the universes are with K = —1, the density matrix again has a strong peak 
at A = 0. The exponent in the density matrix then becomes 

7.3 Summing over topologies 

So far, we have considered the cases that all universes have the same topology. How- 
ever, we can consider a situation where universes with various topologies appear in 
the multiverse. In such case, we should sum over topologies in the multiverse wave 
function. 

To sum over topologies, it is convenient to denote the pair the size and 

the topology of the i-th universe, collectively by Q. Since the probability amplitude 
fi may also depend on the topology of the universe, we denote that with topology ctj 
by fia.., or fiQ. Then, for the multiverse wave function with various topologies, (H2l) is 
generahzed to 

$7v(Ci,--- ,C7v) = JdX (HMo) X <Pe=o{Ci)<Pe=o{C2)---<Pe=o{Cn) w{X)\X). (121) 

We compute the density matrix of our universe from this multiverse wave function. 
By introducing a notation 



dC = y I dz, (122) 



it is given by 



PiC, C) = E / ■ ■ ■ ' (n)^N+i{C, Cl, ■ ■ ■ , ^N) 

= J2m ^(A)Vc'/^c0i^=o(CO*0i^=o(C) X / rfC'Vc"0i^=o(C")l 

OC / (iAw(A)Vc'/^c<^s=o(C')*0E=o(C)exp( E / dz" \iia"(j)E=o{z" ; K, 

J —oo \ „ll J 



(123) 

We note that, compared with the single topology case dS]), the exponent becomes the 
sum over various topologies. 

By comparing (|54ll . (lllSp and (]120p . we find that the flat universes dominate in 
(11231) . Therefore, if universes with any topologies are allowed to emerge, the big fix 
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occurs in such a way that Mmatt in the asymptotic universe with = is minimized. 
In this case the cosmological constant problem is again solved, but the situation for 
the other coupling constants differs much from the case of 5*3 universe. At this stage 
we can not tell which case is more realistic, because we have not specified the details 
about the microscopic dynamics of how universes emerge from nothing with a small 
size z = e. 

8 Summary and Discussion 

In this paper, we have studied the effect of wormholes on the wave function of the 
multiverse and the density matrix of our universe. The wormholes make the multiverse 
wave function a superposition of states with various coupling constants {Aj}. We have 
shown that by examining the density matrix {Aj} are determined in such a way that 
they make the following factor as large as possible 



which is interpreted as the lifetime of the universe. In particular, it is predicted that 
the cosmological constant becomes very close to zero in the far future. If we believe 
the presently observed value of the cosmological constant, which is a non-zero positive 
value, then our analysis suggests that the cosmological "constant" will move towards 
zero such as in the quintessence scenario, where the cosmological constant is the energy 
of a scalar field rolling down in a runaway potential. 

For universes, the coupling constants are determined in such a way that they 
maximize the lifetime of the universe ( |59l) . However, it is difficult to search the max- 
imum point of fl59l) in the parameter space of {Aj} because fl59|) highly depends on 
which parameters are induced by the wormhole effect and also depends on the cosmol- 
ogy and the physics beyond the standard model such as the dark matter and inflation. 
As an illustration of the big fix, we made some assumptions on cosmology and studied 
the possible solution of the gauge hierarchy problem and the strong CP problem. In 
particular, our study suggests that the Higgs mass may be fixed at 



Although we have mainly studied 5*^ universe in this paper, there is a possibility 
that universes are allowed to have the other topologies as in Section [3 We found that 




m 



rrih ~ 



140 ± 20GeV. 



(ma) 
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in such a situation our density matrix is determined only from flat universes, and also 
found that {Aj} are determined such that Mmatt in the far future becomes minimized. 
This naively seems to predict an empty universe and contradict with our universe. 
Therefore, if the universes are allowed to emerge from nothing with any topologies, 
there might be some reason in the quantum gravity that forbids such empty universe 
to emerge as the initial condition. 

In conclusion, the wormhole mechanism is a fascinating scenario since it can solve 
naturalness problems in the standard model and the current cosmology without intro- 
ducing new physics such as supersymmetry or extra dimensions. Although we only 
have presented an illustration of the big flx scenario, it is interesting to explore the pre- 
cise prediction further, and, for this purpose, the deeper understanding of the quantum 
gravity is indispensable. 
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Appendix A Normalization of the Wave Function 

In this appendix, we check the wave function (!37j) satisfles the normalization fl32al) . 



Substituting the wave function, the left hand side is 

/ dz exp(±z / dz\kE>{z')-kE{z'))). (124) 

Jo nz^y kE[zjkE'{z) J 

Note that the delta function can arise from the integral over the asymptotic region 
z ^oo. For large z, kE{z)kE>{z) ~ 9A and kE' - kE ^ fii^' - E) ^ 7fe(^' " 
where we have used A;|; ~ 9A + ^ + ■ ■ ■ . From these, we can check (I124p indeed gives. 




/■°° 1 



exp(±i 



V9A 



1 



{E' - E)\ogz) = 5{E' - E). 
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Appendix B The relation between the curvature 

and e-foldings 



In this appendix, we relate the e-foldings of the initial inflation to the time when 
the curvature term becomes comparable to the energy density. In section |6l we have 
studied the specific case that the curvature term becomes important while the decay 
products of protons are relativistic. We will find that this case corresponds to the 
e-foldings given by f ll3ip . 

We denote by a=„ the scale factor of the universe when ASrad balances with the 
curvature. From (H9|) and is given by 

a. - ASli^^ = {apAMmattf^ (125) 

where AM^att is the total mass of protons in the whole of the universe. It can be 
expressed using the current values of the scale factor qq and the energy density of 
protons Pproton - IGeV/m^: 

AMmatt = CloPproton- (126) 

When the scale factor is around ap, the protons decay, and the decay products, 
especially electrons, are relativistic at that time. However, as the universe expands, 
the energy of these relativistic electrons scales as E electron oc 1/a. And when the scale 
factor becomes about 10^ times as large as a^, they will become non-relativisticc^ 
However, from Conditions in Section [6l the curvature term must become comparable 
to the energy density before it happens. Thus, we have the following constraint on a^, 

ap<a,<apX W\ (127) 

Substituting (11251) and (11261) into the above equation, we obtain 

^ <«o<,/^^^xlOl (128) 



^0 Y Pproton V '^0 V Pproton 

Next, we estimate the ratio ap/a^. Since we have assumed that the cosmological 
constant A is decreasing from the current value to the asymptotic value Kcr — 0, the 
secondary inflation, which is currently going on, ends within a finite time. We denote 



^^The number 10^ comes from a rough estimate of the ratio between the electron mass and its 
energy when it is produced by the proton decay. 
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the e-folding during this inflation by A^. After A gets sufficiently small and the inflation 
ends, the protons dominate the energy density, and the universe scales as a oc 
Thus, Op/ao is given by 

ap/ao~e^lO(=^*^-^°)"i, (129) 

where we have estimated the proton lifetime as Tp ~ lO^^yr and the age of the universe 
today as lO^^yr- Using ffT29D and p^/Z^n - lO^^lyH ffT28D becomes 

^ ^q26/3 < ^0 < ^ ^q26/3+3^ (^3q) 

lO^^ly 

where lO^^ly is the same order as the size of the observable universe and corresponds 
to the lower bound on the e- foldings of the initial inflation, Ne-joid > 55. Thus, the 
above inequality impliej^ 

N N 

- + 75 < N^^foid ^ y + 82. (131) 

Therefore, if Nf._foid is in this range, the cosmological assumption we made in Section 
[6] is realized. 
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